Effects of a Parallel Magnetic Field on the Metal-Insulator Transition in a Dilute 

Two-Dimensional Electron System 
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The temperature dependence of conductivity <r(T) of a two-dimensional electron system in silicon 
has been studied in parallel magnetic fields B. At B = 0, the system displays a metal-insulator 
transition at a critical electron density n c (0), and da/dT > in the metallic phase. At low fields 
(B < 2 T), n c increases as n c (B) — n c (0) oc B 13 (f3 ~ 1), and the zero-temperature conductivity 
scales as a(n s ,B,T — 0)/cr(n s , 0, 0) = f^B 13 /S„) (where 8 n = {n s — n c (0))/n c (0), and n s is electron 
density) as expected for a quantum phase transition. The metallic phase persists in fields of up to 
18 T, consistent with the saturation of n c at high fields. 



PACS Nos. 71.30.+h, 71.27.+a, 73.40.Qv 

The possibility of a metal-insulator transition (MIT) 
in two dimensions (2D) has been a subject of intensive 
research in recent years 0J^] . Most studies have focused 
on the conventional metallic behavior where the conduc- 
tivity a increases with decreasing temperature T (i. e. 
da/dT < 0) but the origin of such a(T) is still not under- 
stood. Similarly, a dramatic decrease of a with magnetic 
fields B applied parallel to the 2D plane remains puz- 
zling. Above some characteristic field, which is a function 
of the carrier density n s , this unexpectedly large negative 
magnetoconductance is followed by a weaker dependence 
on B. Parallel B has been also shown to suppress the 
metallic T-dependence (da/dT < 0) and transform 
it into an insulatinglike behavior [da/dT > 0). Since it 
has been commonly assumed that da/dT > is necessar- 
ily related to an insulating state at T = 0, this was taken 
as evidence that high parallel B eliminates the MIT and 
the 2D metallic phase entirely OJ . The field where da / dT 
changes sign (for a fixed n s ) has been identified as the 
critical field for a field- induced MIT j5j. Furthermore, it 
has been established that this critical field is compara- 
ble both to the field where the crossover from a low-field 
to high-field regime takes place and to the field where 
2D carriers become fully spin-polarized suggest- 
ing that the 2D metallic phase may exist only in a spin- 
unpolarized system. By contrast, the methods that do 
not use the sign change of da/dT as a criterion to deter- 
mine n c have yielded || finite values of n c even in high 
parallel B where da / dT > at all n s , which implies that 
the 2D metal may exist even for spinless electrons. It is 
clear that the fate of the metallic phase in a parallel B 
represents one of the major open issues in the studies of 
dilute, strongly interacting 2D systems. 

The main difficulty associated with the attempts to re- 
solve this issue so far has been that a(T) does not 
have a simple form for n s > n c (0) in parallel B. There- 
fore, this does not allow one to make reliable extrapola- 
tions to T — and to establish the existence of a true 



MIT with high credibility. In addition, there is no theo- 
retical justification for assuming that the 2D metallic 
phase can be characterized only by da/dT < 0. In dis- 
ordered 3D metals, for example, it is well known that 
this derivative can be either negative or positive near the 
MIT fll|| . More relevant to this work, a recent study j^j 
of a 2D electron system in Si metal-oxide-semiconductor 
field-effect transistors (MOSFETs) in B = has estab- 
lished unambiguously that the metallic phase may be de- 
scribed by da/dT > such that a decreases to a nonzero 
value as T — > 0. This occurs when some of the localized 
states in the tail of the upper electric subband |l2j also 
become populated 14 15|] as the subband splitting is re- 
duced by applying voltage (V su b) to the Si substrate [ p^[ . 
Such localized states act as additional scattering centers 
for 2D electrons p6| , and perhaps may even act as local 
magnetic moments if they are singly occupied (due to a 
large on-site Coulomb repulsion) |l5|Jl7| ]. 

Regardless of the detailed microscopic picture, the 2D 
metal with da/dT > in B — exhibits two fea- 
tures that distinguish it from the da/dT < case, 
and that are of crucial importance for the study of a 
in parallel B. First, a(T) follows a simple and pre- 
cise form over a very broad (two decades) range of T: 
a(n s ,T) = a(n s ,T = 0)+C(n s )T 2 . This allows a reliable 
extrapolation to T = 0, which yields a(n s ,T = ) ~ 8fi 
(fj, w 3), as expected in the vicinity of the MIT 0!§. 
Second, near the MIT, the data obey dynamical scaling 
a{n s ,T) = a c {T)f(T/d^) with a c = a{n s = n c , T) ~ T x 
(zv = 1.3 ± 0.1, x w 2.6, fi = x(zv) = 3.4 ± 0.4), 
both in agreement with theoretical expectations near a 
quantum phase transition JTo| | and consistent with the 
extrapolations of a(T) to T — 0. Here we show that 
such features allow us to establish a phase diagram in 
the (6 n ,B,T = 0) plane. We find that the behavior of 
low-field a(B) is consistent with the phase diagram, and 
can be attributed to an increase in n c (B), as expected 
near a true MIT p0[. The saturation of n c (B) observed 
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FIG. 1. o(T) for sample 9 at B — 4 T. n a varies from 
3.0 x W rl cm~ 2 (top) to 0.7 x W rl cm' 2 (bottom) in steps 
of 0.1 x 10 n cm- 2 . n c (B = 4 T) = 1.4 X lO^cm -2 and 
is marked by the arrow. The dashed lines guide the eye. <r c 
clearly follows a simple power-law dependence on T : o c oc T x . 



at higher B strongly suggests that the 2D metal may 
exist even for spinless electrons. 

Measurements were carried out on n-channel Corbino 
shaped Si MOSFETs (channel length = 0.4 mm, 
mean circumference = 8 mm) with poly-Si gates, self- 
aligned ion- implanted contacts, and the peak mobility 
~ 1 m 2 /Vs at 4.2 K. Other sample details have been 
given elsewhere Jj^JToj pl . n s was controlled by the gate 
voltage V g . V su b = +1 V was applied to minimize 
the subband splitting and maximize the T-range (up to 
~ 2 K) in which the metallic behavior with da/dT > is 
observed fL5|,3 ■ Conductance was measured as a function 
of V g using standard low-noise analog lock-in techniques 
at ^ 17 Hz and a low- noise current preamplifier (the lead 
resistance was subtracted by the usual method). Exci- 
tation voltages were kept low enough to ensure that the 
conduction was Ohmic at the lowest T (see Ref. [2] for 
more details). A systematic study of a(n s ,T) was done 
in a He 3 cryostat in static parallel fields up to 9 T and for 
0.25 < T < 4 K. Measurements up to 18 T were carried 
out with the sample in the mixing chamber of a dilution 
refrigerator with a base T « 0.020 K. 

Figure [l] shows cr(T) for sample 9 (also studied in 
Ref. [2]) at B = 4 T for various n s in the vicinity of 
the MIT. The data are qualitatively similar to the .8 = 
case Q]. For example, at the lowest n s (i. e. n s < n c ), 
a decreases exponentially with decreasing T, indicating 
an insulating state at T = 0. In the metallic phase 
(n s > n c ), a(T) is weaker and its curvature is the op- 
posite from the one expected for an insulating state. It 
clearly extrapolates to a finite value as T — > 0, as dis- 
cussed in more detail below. As shown in Fig. [I], n c is 
identified as the density where a c — a(n s — n c , T) oc T x , 
consistent with the B = case and in agreement with 
general arguments po| . The exponent x is found to be 
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FIG. 2. T — phase diagram for two samples. The dashed 
line guides the eye. The boundary between metallic and insu- 
lating phases is described by a power-law relation (see inset) 
[n c {B) -n c (0)]/ra c (0) oc (B/B*) at low fields, with the same 
crossover exponent j3 ~ 0.9 for both samples (B* = 1 T for 
sample 9). Inset: the same data vs. B on a log-log scale. The 
dashed lines are fits with the slopes equal to (3. At B — 0, 
n c (10 n cm- 2 ) = 0.95±0.05 and 0.85±0.05 for samples 9 and 
17, respectively. 

2.7 ± 0.4, and remains constant as a function of B. 

The critical density n c was determined in this way at 
each given B. Fig. || inset shows the relative change 
of n c (B) with respect to its zero-field value n c (Q) as a 
function of B for two different samples. At low fields 
(B < 2 T), n c (B) increases with B in a power-law fash- 
ion: [n c (B) - n c (0)]/n c (0) oc B@, where the crossover 
exponent (5 = 0.80 ± 0.06 and 0.9 ± 0.1 for samples 9 
and 17, respectively. We point out that such a power- 
law shift of n c with B is expected to occur in the case 
of a true MIT p0| , and has been observed in several 3D 
systems In the case of conventional 2D metallic 

behavior in high-mobility Si MOSFETs, finite values of 
n c (B) were found based on vanishing activation energy 
and vanishing nonlincarity of current-voltage curves as 
extrapolated from the insulating side |J . For those data, 
we find that the corresponding f3 = 1.1 ± 0.1. It is in- 
teresting that the crossover exponents /3 are almost the 
same in both cases (Ref. || and Fig. ||) even though the 
B = behaviors of da/dT are very different. Further- 
more, in both cases the saturation of n c (B) is observed 
at higher fields. The determination of n c (B) based on a 
sign change of da/dT has yielded (3 « 1 in high- mobility 
Si MOSFETs [|, and (3 w 1.4 in a 2D hole system in 
GaAs H . While the difference in (3 between Si and GaAs 
devices might be due to the difference in the dominant 
scattering mechanisms in the two materials, it should be 
noted that the mere sign of da/dT is clearly not a reliable 
indicator of the nature of the (T = 0) ground state. 

Since (3 ~ 0.9 for both of our samples and only the 
absolute values of n c (B) are different, it is possible to 
create a generalized phase diagram in the (5 n ,B,T — 



2 



0) plane as shown in Fig. |[ The dashed line in Fig. || 
represents the boundary between metallic and insulating 
phases. It is clear that, for a fixed, small S n , the system 
will undergo a magnetic field-driven MIT, leading to a 
strong suppression of magnetoconductance (MC). On the 
other hand, for densities far from n c (0) (i. e. large S n ), 
the 2D system will remain metallic with increasing B. 
In that case, one expects a much weaker MC. Exactly 
this kind of behavior of MC has been observed in both 
our samples [|o) and high-mobility Si MOSFETs in the 
conventional metallic regime S. 

First we focus on the regime of low fields (B < 2 T) and 
small <5„, where MC exhibits a large drop |^(|. We find 
that all the MC curves a(B)/a(0) for a fixed n s and dif- 
ferent T can be collapsed onto one function using a single 
scaling parameter T(n s , T) [Fig.||(a)]. For each given n s , 
we find empirically that T{n s ,T) = a{n s ) + b{n s )T 1 - z±0 - 2 
[Fig. 0(b)]. The parameter a(n s ) decreases with decreas- 
ing n s and extrapolates to zero at a density equal to 
(0.88 ± 0.04) x 10 n cm~ 2 . Within our measurement er- 
ror, this density is the same as n c (0). Fig. ||(a) inset 
shows that a fit to the critical form, a(n s ) = ao5l/^, 
yields a crossover exponent (3 = 1.2±0.2 (ao = 1.7±0.4), 
consistent with the value of (3 that was established inde- 
pendently earlier (Fig. ||) based on the extrapolations 
of cr(T) to T = in the metallic regime at fixed B 
(Fig. [I]). Moreover, Fig. ||(a) shows that it is also possi- 
ble to collapse a(B)/a(0) for all different n s and T using 
a single scaling parameter Bq(ti s ,T) = A(n s )J r (n s ,T), 
which combines the previous scaling parameter T with 
a density dependent prefactor A(n s ). The fitting pa- 
rameters A(n s ) = (2.8 ± 0.3)n s - (1.6 ± 0.3) and b(n s ) = 
-(11.4±0.9)n s + (14.9±0.9) (all n s in units of 10 n cm- 2 ) 
do not show much variation near n c (0). 

Therefore, our results demonstrate that the zero- 
temperature conductivity scales with n s and B as 



c(<5„,-B,0)/ct((5„,0,0) = f(B/5l/P), 



(1) 



where (3^1. For a fixed n s , the crossover function 
/(B/Sn 13 ) decreases by orders of magnitude with increas- 
ing B [Fig. ||(a)] as the system approaches the magnetic 
field-driven MIT. The critical fields B c for such a tran- 
sition are given by B c oc 5n in agreement with the 
phase diagram determined earlier (Fig. [2]). It is clear 
that the scaling parameter Bq(ti Si T — 0) oc B c . We 
also note that the scaling plot in Fig. ||(a) includes a 
few densities below n c (0) in the quantum critical re- 
gion, where T{n s ,T) oc T 13 in the given T range. 
Furthermore, for n s — n c (0), our results show that 
a(0,B,T)/a(0,0,T) = fjB/T 1 - 3 ). According to stan- 
dard scaling arguments |l0[ , the exponent 1.3 = (3/zv. 
With f3 = 1.2 ± 0.2, this yields zv = 0.9 ± 0.3, in rea- 
sonable agreement with zv — 1.3 ± 0.1 determined from 
scaling of a(n s ,T) in B — All these findings pro- 
vide strong experimental evidence for a quantum phase 
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FIG. 3. (a) Scaling of low-field (B < 2 T) magne- 
toconductivity for n 3 (10 11 cm -2 ) = 0.80 — 1.10 in steps 
of 0.02, and ten different T from 0.25 to 1.4 K. The 
scaling parameter Bo(n s ,T) = A(n a )T(n a , T). (b) 
T(n s ,T) — a(n s ) + 6(n s )T 1 ' 3 , as shown for several n s . 
a(n s ) oc ' , as shown in the inset of (a). 



transition in both B — and B ^ in our system. 

Recent studies of some Si MOSFETs in the conven- 
tional metallic regime have shown |2lJ that, in that case, 
the MC curves can be also scaled with a single parameter 
B a , where B a corresponds to the onset of full spin po- 
larization of the electron system |22j {[3 « 1.7 and (3 w 1 
in Refs. ||l[] and [^2), respectively). However, since the 
form of o~(T) in the conventional metallic regime does 
not allow reliable extrapolations to T = in B ^ 0, it 
is difficult to establish the existence of a metallic phase 
and n c (B) (or, equivalently, B c (n s )) with certainty. We 
showed earlier @,|l§ that our Si MOSFETs exhibit ex- 
actly the same kind of conventional metallic behavior in 
absence of scattering by local magnetic moments. There- 
fore, we expect that similar values of B will be required 
for full spin polarization of our 2D system at compara- 
ble n s , even when a small fraction of electrons occupies 
the tail of the upper subband. Our work demonstrates 
that, in spite of the qualitative differences in o~(T) in 
the two regimes (with and without local magnetic mo- 
ments), the effect of the parallel B on a is remarkably 
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FIG. 4. Sample 19, B = 18 T: <r(T) plotted vs. T 3 for 
n s (10 11 cm -2 ) from 4.9 (top) to 3.8 (bottom) in steps of 0.1; 
0.020 < T < 0.478 K. The dashed lines are fits a fit (n a ,T). 
Inset: The deviation of the data from the fit vs. temperature. 



similar. This strongly suggests that in our system at 
low fields (B < 2 T) , the critical fields B c ~ B a , where 
B a is the polarization field. Indeed, in that regime B c 
in Fig. H is comparable to B a in Ref. |^] for the same 
n s . In addition, since B a oc 5 n continues to hold over a 
wide range of n s (up to 5 n ~ 10) P4] , this implies that 
for <5 n ~ (0.5 — 0.6) in Fig. |2|, the 2D electrons will be 
fully spin polarized at ~ 2 T. Therefore, at B > 2 T 
the MIT will occur between a spin-polarized metal and 
a spin-polarized insulator. This conclusion is supported 
by our observation that n c no longer seems to depend on 
the magnetic field for B > 2 T. 

The detailed properties of the apparent spin-polarized 
2D metallic phase and the concomitant MIT require fur- 
ther careful investigation that is beyond the scope of this 
Letter. Here we demonstrate only that the metallic phase 
persists even in fields as high as 18 T, where we esti- 
mate [^2) that the 2D system is fully spin polarized in 
the range of n s studied. We find that a(T) in the metal- 
lic phase follows a simple power-law form 



a(n s , B, T) = a{n s ,B, T = 0) + C(n s , B) T 



'(B) 



(2) 



similar to the B — case, where a(B = 0) = 2.0 ±0.1 [|. 
In both cases, finite values of a(n s ,B,T = 0) indicate 
that the 2D system is in the metallic state. At B = 18 T, 
however, a(T) is best described with a — 3 in Eq. (j^), 
as shown in Fig. ||. While the mechanism that leads to 
a = 3 is unclear, it is interesting that the T 3 term (albeit 
with a different sign) has been observed in the quasi- 
2D metallic Sr3Ru207 below 0.5 K near a metamagnetic 
quantum critical point p3[ . In our case, such a jump 
in magnetization would occur as 2D electrons rearrange 
from an unpolarized to a spin-polarized metallic state 
with decreasing n s at a fixed high B. The T 3 variation of 
resistivity has been also observed in a low carrier density 
metallic state of SrB 6 p4[ but the microscopic mechanism 



is still not understood. 

In summary, we have studied the effect of a paral- 
lel B on the conductivity of a 2D electron system with 
da/dT > in the metallic phase, even at B = 0. The 
zero-temperature (8 n ,B) phase diagram has been con- 
structed based on the simple form of cr(T) in the metal- 
lic regime. Furthermore, the low-field magnetoconduc- 
tivity has been found to scale with n s , consistent with 
the (8 n ,B,T = 0) phase diagram. These findings pro- 
vide strong experimental evidence for a quantum phase 
transition in both B = and B ^ 0. The metallic phase 
is found to persist in fields as high as 18 T. Our results 
strongly suggest that the 2D metallic phase may exist 
even in the absence of spin degrees of freedom. 
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